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We suggest that tunable orientational nonlinearity of nematic liquid crystals can 
be employed for all-optical switching in periodic photonic structures with liquid- 
crystal defects. We consider a one-dimensional periodic structure of Si layers with a 
local defect created by infiltrating a liquid crystal into a pore, and demonstrate, by 
solving numerically a system of coupled nonlinear equations for the nematic director 
and the propagating electric field, that the light-induced Freedericksz transition can 
lead to a sharp switching and diode operation in the photonic devices. 

I. INTRODUCTION 

During the past decade, photonic crystals (artificially fabricated one-, two- and three- 
dimensional periodic dielectric materials) have attracted a great deal of interest due to 
their ability to inhibit the propagation of light over special regions known as photonic band 
gaps [1]. Such photonic bandgap materials are expected to revolutionize integrated optics 
and micro-photonics due to an efficient control of the electromagnetic radiation they provide, 
in such a way as semiconductors control the behavior of the electrons Q|. 

In general, the transmission of light through photonic crystals depends on the geometry 
and the index of refraction of the dielectric material. Tunability of the photonic bandgap 
structures is a key feature required for the dynamical control of light transmission and various 
realistic applications of the photonic crystal devices. One of the most attractive and practical 
schemes for tuning the band gap in photonic crystals was proposed by Busch and John (^J, 
who suggested that coating the surface of an inverse opal structure with a liquid crystal could 



be used to continuously tune the band gap, as was confirmed later in experiment 



. This 
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original concept generated a stream of interesting suggestions for tunable photonic devices 
based on the use of liquid crystals infiltrated into the pores of a bandgap structure 0. 
The main idea behind all those studies is the ability to continuously tune the bandgap 
spectrum of a periodic dielectric structure using the temperature dependent refractive index 
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of a liquid crystal |4|, y, |y, ML or its property to change the refractive index under the action 
of an applied electric field ja 0, 13 • 

Another idea of the use of liquid crystals for tunability of photonic crystals is based 
on infiltration of individual pores and creation of liquid crystal defects Il3, Q> 



and even defect-induced waveguide circuits |1 1| . In this case, the transmission properties 
can be controlled, for example, by tuning resonant reflections associated with the Fano 
resonances Q] observed when the frequency of the incoming wave coincides with the 
frequency of the defect mode. As a result, the defect mode becomes highly excited at the 
frequency of the resonant reflection, and it can be tuned externally, again by an electric field 
or temperature. 

However, liquid crystals by themselves demonstrate a rich variety of nonlinear phenomena 
(see, for example, Refs. 12, 1£, 3,13)- Therefore, nonlinear response of liquid crystals can 
be employed for all-optical control of light propagation in periodic structures and tunability 
of photonic crystals. In this paper, for the first time to our knowledge, we analyze the 
possibility of tunable all-optical switching in one-dimensional periodic structure with a liquid 
crystal defect. We demonstrate that the light field with the intensity above a certain critical 
value corresponding to the optical Freedericksz transition changes the optical properties 
of the liquid- crystal defect such that the nonlinear transmission of the photonic structure 
allows for all-optical switching, and the similar concept can be employed for creating of a 
tunable all-optical diode. 



II. NONLINEAR TRANSMISSION OF A LIQUID CRYSTAL SLAB 

First, we study the light transmission of a single slab of nematic liquid crystal and derive 
a system of coupled nonlinear equations for the liquid-crystal director reorientation in the 
presence of the propagating electric field of a finite amplitude. The corresponding steady- 
state equation for the director n can be obtained by minimizing the free energy which can 



be written in the following form 13, |21 ] 
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where / e i is the elastic part and f opt is the optical part of the energy density. Here #n, #22 
and #33 are splay twist, and bend elastic constants, respectively, D = eE, e is the dielec- 
tric tensor, and the real electric field is taken in the form E rea ; = (l/2)[E(r) exp(— iut) + 
E*(r)exp(iwt)]. 

We assume that linearly polarized light wave propagates normally to the liquid-crystal 
slab with the initial homeotropic director orientation along z [see Fig. IHa)]. Under the 
action of the electric field polarized outside the slab along x, the director can change its 
direction in the (x, z) plane and, therefore, we write the vector components of the director 
in the form n = {sin0(z), 0, cos0(,z)}. Then the elastic part of the free energy density can 
be written as 

\ ■ (2) 
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Taking into account that the dielectric tensor e can be expressed in terms of the director 
components, = e±5ij + e a ninj, where e a = e\\— e± and e||, e± are the liquid crystal dielectric 
constants at the director parallel and perpendicular to the electric vector, we can write 

e a sin cos <p \ 
e± 

V e a sin cos e± + e a cos 2 J 

As a result, the optical part of the free energy density takes the form 



^ e_L + e a sin 2 . 
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After minimizing the free energy (JIJ) with respect to the director angle 0, we obtain the 
nonlinear equation for the director in the presence of the light field 
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FIG. 1: Nonlinear transmission of a liquid-crystal slab, (a) Schematic of the problem. (b,c) 
Maximum angle of the director and transmission vs. the light intensity in the slab. Blue and red 
curves correspond to the increasing and decreasing light intensity, respectively. 

where A(<f>) = (K n sm 2 (p + K 33 cos 2 0), B((p) = (K 33 — Ku) sin cos 0, and we take into 
account that, as follows from D z = 0, that the electric vector of the light field has the 
longitudinal component, E z = — (e xz /e zz )E x = — [e a sin0cos0/(e_i_ + e a cos 2 (p)]E x (see also 
Ref. Q). 

From the Maxwell's equations, we obtain the equation for the electric field E x , 

d 2 Er 



dz 2 



e ± (e± + e a ) 
e ± + e a cos 2 



-E x = 0, 



(5) 



where k = 2n/ (Ac). Moreover, it can be shown that the z component of the Poynting 

vector I = S z = (c/8n)E x H* remains constant during the light scattering and, therefore, it 
can be used to characterize the nonlinear transmission results. 

As the boundary conditions for the coupled nonlinear equations (J1J) and ©, we assume 
that there is an infinitely rigid director anchoring at both surfaces of the slab, i.e. 



0(0) = <j>{L) = 0, 



(6) 



and also introduce the scattering amplitudes for the optical field 



EJz) 



S^exp^ikz) + £ re {exp(—ikz), z < 0, 
S OVLt exp(ikz), z>L, 



(7) 



where L is the thickness of the liquid-crystal slab, £ in , £ re f, and £ out are the electric field 
amplitudes of incident, reflected, and outgoing waves, respectively. 

To solve this nonlinear problem, first we fix the amplitude of the outgoing wave S out and 
find unique values for the amplitudes of the incident, & mi and reflected , £ re f, waves. By 
using the so-called shooting method [3], in Eq. (p|J) for the director we fix the amplitude of 
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the outgoing wave and, assuming that <f>(L) = at the right boundary, find the derivative 
(d<f)/dz) z= L such that after integrating we obtain a vanishing value of the director at the 
left boundary, i.e. 0(0) = 0. Because Eq. (j3J) is a general type of the nonlinear pendulum 
equation, we look for periodic solutions with the period 2L. Obviously, there exists an 
infinite number of such solutions and, therefore, there is an infinite set of the derivatives 
(d(p/dz) z= L which satisfy Eq. (J3J) and the condition JBJ). All such solutions correspond to 
some extrema points of the free energy of the system. However, we are interested only 
in that solution which realizes the minimum of the free energy. By analyzing our coupled 
nonlinear equations in a two-dimensional phase space, we can show that the corresponding 
solution lies just below the separatrix curve, and it has no node between the points z = 
and z = L. This observation allows us to reduce significantly the domain for our search for 
the required values of the derivative (dcf>/dz) z= L- 
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FIG. 2: Transmission of an one-dimensional periodic structure with an embedded liquid-crystal 
defect. In the linear regime, the transmission is characterized by the presence of an in-gap resonant 
peak due to the excitation of a defect mode. Nonlinear transmission displays bistability at the 
defect-mode frequency with two different thresholds for " up" and " down" directions and a hysteresis 
loop (see the insert). 



The obtained solutions can be characterized by the maximum angle max of the director 
deviation which, as is intuitively clear, should be reached near or at the middle of the slab. 
In Fig.^b,c), we plot the maximum angle max and the transmission coefficient of the liquid- 
crystal slab, defined as T = |£ ut| 2 /|£in| 2 , vs. the light intensity. For numerical calculations, 
we use the parameters K x \ = 4.5 x 10~ 7 dyn, K 33 = 9.5 x 10~ 7 dyn, e a = 0.896, e± = 2.45, 



L = 200nm, and A = 1.5/im, that correspond to the PAA liquid crystal 



231 ]: because of a 
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lack of the corresponding data at the wavelength A = 1.5/im, the values of the dielectric 
constant are taken from the optical range. 

From the results presented in Fig.EJb,c), we observe sharp jumps of the director maximum 
angle max and the transmission coefficient T due to the optical Freedericksz transition in 
the liquid-crystal defect. However, a variation of the transmission coefficient during this 
process is not larger than 15%. The threshold of the optical Freedericksz transition appears 
to be different for the increasing and decreasing intensity of the incoming light, so that this 
nonlinear system is bistable, and it displays a hysteresis behavior. The bistable transmission 
of the liquid-crystal slab is similar to that predicted for the slab of PAA liquid crystal in 
the geometric optics approximation Q], and such a behavior is explained by the existence 



of the metastable state which the system occupies at the decreasing light intensity [Tj, lfll ]. 



III. LIQUID-CRYSTAL DEFECT IN A PERIODIC PHOTONIC STRUCTURE 

Now, we study the similar problem for a liquid-crystal defect infiltrated into a pore of 
the periodic structure created by Si layers with the refractive index n = 3.4. For simplicity, 
we consider a one-dimensional structure with the period a = AOOnm and the layer thickness 
d\ = 200nm, that possesses a frequency gap between \A\im and 2.5/zm. We assume that one 
of the holes is infiltrated with a PAA nematic liquid crystal with e± = 2.45. Such a defect 
modifies the linear transmission of the periodic structure by creating a sharp defect-mode 
peak at the wavelength A^ ~ 1.5fim, as shown in Fig. El 



To solve the nonlinear transmission problem, we employ the transfer matrix approach [24 1 
implementing it for the solution of the full system of coupled equations (j3J) and ©. By 
tuning the input intensity at the defect mode, we observe the same scenario as for a single 
liquid-crystal slab [cf. the insert in Fig. Eland Fig. ^c)]. Namely, there exists a hysteresis 
loop in the transmission with two different thresholds for the increasing and decreasing 
intensities. The difference is, however, in the values. Due to a small width of the defect- 
mode resonance, even a small reorientation of the director leads to a sharp (up to 90%) 
change in the transmission. Another significant difference is that the threshold values are 
lower by four orders of the magnitude, for a given periodic structure for which we take 10 
layers from each side of the defect. 

Finally, we notice that in a finite periodic structure the defect placed asymmetrically 
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FIG. 3: Example of a tunable all-optical diode based on the optical Freedericksz transition in a 
liquid-crystal defect. Asymmetrically placed defect leads to different threshold intensities of the 
switching for the waves propagating from the right and left, respectively. 

(see Fig. EJ) allows to create a nonreciprocal device when the threshold intensities for the 
molecular reorientation differ for the light propagating from the right and left. This feature 
is associated with the operation of an optical diode HQ- As can be seen in Fig. El by 
shifting the infiltrated liquid-crystal defect closer to one of the edges of the structure and 
fixing the total length of the structure, we can increase the switching power and extend the 
diode operation region decreasing the transmission power. Moreover, these results show that 
the threshold intensities depend strongly on the number of periods to the structure edge, 
due to a stronger confinement of the defect mode. Also, this gives us a possibility to reduce 
significantly the switching power simply by taken larger number of periods in the photonic 
structure. 
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IV. CONCLUSIONS 

We have demonstrated that the orientational nonlinearity of nematic liquid crystals can 
be employed to achieve tunable all-optical switching and diode operation in periodic photonic 
structures with infiltrated liquid-crystal defects. For the first time to our knowledge, we have 
solved a coupled system of nonlinear equations for the nematic director and the propagating 
electric field for the model of a one-dimensional periodic structure created by Si layers with a 
single (symmetric or asymmetric) pore infiltrated by a liquid crystal. We have demonstrated 
that the threshold of the optical Freedericksz transition in the liquid-crystal defect is reduced 
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dramatically due to multiple reflections in the periodic structure, so that such a defect may 
allow a tunable switching and diode operation in the photonic structure. 

Acknowledgements 

The work has been supported by the Australian Research Council. The authors thank 
B.Ya. Zeldovich, I.C. Khoo, M. Karpierz, and O. Lavrentovich for useful discussion of our 
results and suggestions, and I.V. Shadrivov for the help in numerical simulations. 



[1] J.D. Joannopoulos, R.D. Meade, and J.N. Winn, Photonic Crystals: Molding the Flow of 

Light (Princeton University Press, Princeton, NY, 1995). 
[2] E. Yablonovitch, "Inhibited spontaneous emission in solid-state physics and electronics," Phys. 

Rev. Lett. 58, 2059-2062 (1987). 
[3] K. Busch and S. John, "Liquid-crystal photonic-band- gap materials: The tunable electromag- 
netic vacuum," Phys. Rev. Lett. 83, 967-970 (1999). 
[4] K. Yoshino, Y. Shimoda, Y. Kawagishi, K. Nakayama, and M. Ozaki, "Temperature tuning 

of the stop band in transmission spectra of liquid-crystal infiltrated synthetic opal as tunable 

photonic crystal," Appl. Phys. Lett. 75, 932-934 (1999). 
[5] S.W. Leonard, J.P Mondia, H.M. van Driel, O. Toader, S. John, K. Busch, A. Birner, U. 

Gosele, and V. Lehmann, "Tunable two-dimensional photonic crystals using liquid-crystal in- 
filtration," Phys. Rev. B 61, R2389-R2392 (2000). 
[6] Ch. Schuller, F. Klopf, J.P. Reithmaier, M. Kamp, and A. Forchel, "Tunable photonic crystals 

fabricated in III-IV semiconductor slab wavelengths using infiltrated liquid crystals," Appl. 

Phys. Lett. 82, 2767-2769 (2003). 
[7] T. T. Larsen, A. Bjarklev, D. S. Hermann, and J. Broeng, "Optical devices 

based on liquid crystal photonic bandgap fibers," Opt. Express 11, 2589-2596 (2003), 

|http://www.opticsexpress.org/abstract.cfm?URI=OPEX-ll-20-2589 
[8] D. Kang, J.E. Maclennan, NA. Clark, AA. Zakhidov, and R.H. Baughman, "Electro-optic 

behavior of liquid-crystal-filled silica photonic crystals: Effect of liquid-crystal alignment," 

Phys. Rev. Lett. 86, 4052-4055 (2001). 



9 



[9] M.J. Escuti, J. Qi, and G.P. Crawford, "Tunable face-centered-cubic photonic crystal formed 
in holographic polymer dispersed liquid crystals," Opt. Lett. 28, 522-524 (2003). 
[10] E. Graugnard, J.S. King, S. Jain, C.J. Summers, Y. Zhang- Williams, and I.C. Khoo, "Electric- 
field tuning of the Bragg peak in large-pore TiO"2 inverse shell opals," Phys. Rev. B 72, 
233105-4 (2005). 

[11] S.F. Mingaleev, M. Schillinger, D. Hermann, and K. Busch, "Tunable photonic crystal circuits: 
concepts and designs based on single-pore infiltration," Opt. Lett. 29, 2858-22860 (2004). 

[12] I. Del Villar, I. R. Matias, F. J. Arregui, and R. O. Claus, "Analysis of one- 
dimensional photonic band gap structures with a liquid crystal defect towards devel- 
opment of fiber-optic tunable wavelength filters," Opt. Express 11, 430-436 (2003), 
http://www.opticsexpress.org/abstract. cfm?URI=OPEX-ll-5-430 

[13] R. Ozaki, T. Matsui, M. Ozaki, and K. Yoshino, "Electrically color-tunable defect mode lasing 
in one-dimensional photonic band-gap system containing liquid crystal," Appl. Phys. Lett. 82, 
3593-3595 (2003). 

[14] E.P. Kosmidou, E.E. Kriezis, and T.D. Tsiboukis, "Analysis of tunable photonic crystal de- 
vices comprising liquid crystal materials as defects," IEEE J. Quantum Electron. 41, 657-665 
(2005). 

[15] S. Fan, "Sharp asymmetric line shapes in side-coupled waveguide-cavity systems," App. Phys. 

Lett. 80, 908-910 (2002). 
[16] A.E. Miroshnichenko and Y.S. Kivshar, "Sharp bends in photonic crystal 

waveguides as nonlinear Fano resonators," Opt. Express 13, 3969-3976 (2005), 

http:/ /www.opticsexpress.org/abstract.cfm?URI=OPEX-13-ll-3969 
[17] B.Ya. Zel'dovich, N.V. Tabiryan, and Yu.S. Chilingaryan, "Freedericksz transition induced by 

light fields," Zh. Eksp. Teor. Fiz. 81, 72 (1981) [Sov. Phys.-JETP 81, 72 (1981)]. 
[18] I.C. Khoo, "Optically induced molecular reorientation and third order nonlinear processes in 

nematic liquid crystals," Phys. Rev. A 23, 2077-2081 (1981). 
[19] H.L. Ong, "Optically induced Freedericksz transition and bistability in a nematic liquid crys- 
tal," Phys. Rev. A 28, 2393-2407 (1983). 
[20] N.V. Tabiryan, A.V. Sukhov, and B.Ya. Zel'dovich, "Orientational optical nonlinearity of 

liquid crystals," Mol. Cryst. Liq. Cryst. 136, 1-131 (1986). 
[21] P.G. de Gennes, The Physics of Liquid Crystals, (Clarendon Press, Oxford, 1979). 



10 



[22] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes in C++, 

(Cambridge University Press, 2002). 
[23] M.J. Stephen, and J.P. Straley, "Physics of Liquid Crystals," Rev. Mod. Phys. 46, 617-704 

(1974). 

[24] P. Yeh, Optical Waves in Layered Media, (John Wiley & Sons, New York, 1988). 

[25] M. Scalora, J.P. Dowling, CM. Bowden, and M.J. Bloemer , "The photonic band edge optical 

diode," J. App. Phys. 76, 2023-2026 (1994). 
[26] K. Gallo, G. Assanto, K.R. Parameswaran, and M.M. Fejer, "All-optical diode in a periodically 

poled lithium niobate waveguide," App. Phys. Lett. 79, 314-316 (2001). 



